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Abstract. Muon cooling for a neutrino factory or muon collider can be achieved using 
low-Z absorbers in strong focusing fields. Proposed cooling lattices place absorbers in 
solenoidal fields ranging up to 30 to 40 T. The cooling performance of these lattices is 
determined by the interplay of ionization energy loss and Moliere scattering, but Bet he's 
classic treatment of Moliere scattering ignores the helical motion of charged particles in 
solenoidal fields. When this motion is taken into account, the performance of these lattices 
can be better than predicted by simulations using the standard treatment. 

Contribution to NUFACT 11, Xlllth International Workshop on Neutrino Factories, Super 
beams and Beta beams, 1-6 August 2011, CERN and University of Geneva 
(Submitted to lOP conference series) 



1. Introduction 

The question of the effect of strong magnetic fields on multiple scattering was raised by 
Lebrun in 1999 [Ij. He realized that Bethe's derivation [2J of the Moliere scattering 
distribution implicitly assumed straight-line motion between scatters. In contrast, in a 
strong magnetic field, a charged particle executes circular or helical motion, and in the 
limit of infinitely strong field, such a particle's path is confined to a magnetic-field line. 
Clearly, this must reduce the multiple-scattering-induced spread in position as a charged 
particle traverses a scattering medium, as Lebrun demonstrated using Geant 3. Since this is 
exactly the configuration proposed in many implementations of muon cooling, simulations 
of muon-cooling channels risk underestimating their effectiveness, especially as the beam 
emittance shrinks and the field strength grows to 20 T and beyond [3J. 

In order to understand the degree of underestimation that may occur in such simulations, 
we undertook to compare the treatment of scattering in magnetic fields in the principal 
codes used in recent MAP and NFMCC cooling studies: G4beamline [4J and ICOOL [5j. 
To provide a simple and well-understood baseline, we also evaluated the effect using a 
Gaussian approximation to the scattering distribution]^ 

2. Gaussian Model 

In our Gaussian study, we let each muon start at one end {z = 0) of a 30-cm-long liquid- 
hydrogen scatterer (a typical muon-cooling-channel absorber length), with longitudinal 



^ For simpUcity, we disabled dE/dx energy loss in all three calculations. 
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Figure 1. Distributions in (left) x and (right) tan^^^ at the end of a 30 cm LH2 absorber 
in Gaussian multiple-scattering model, for (top) 5 and (bottom) 100 absorber slices. 



momentum pz — 200MeV/c and transverse momenta Px — Py — 0. We subdivided the 
absorber into thinner and thinner longitudinal slices, tracing each muon's path as a straight 
line (if B — 0) oy helix (if S 7^ 0) between the centers of successive slices. At the center of 
each slice, we allowed 9x and 9y to change by random amounts, with r.m.s. variation [6j 

(1) 

where Az is the width of a slice and Xq — 866 cm is the radiation length of liquid 
hydrogen [6J. Figure [l] shows an example of the transverse spatial and angular distributions 
at the end of the absorber for the cases 5 = and ngiices = 5 or 100. As expected, the 
number of slices makes no difference if S = 0. 

We then allowed the longitudinal magnetic field to increase up to 100 T and studied 
the number of slices required for the r.m.s. x value at the end of the absorber to stabilize, 
finding that 10 slices suffice (Fig. [2]). We observe a dramatic suppression of spatial (but not 
angular) spread as the magnetic field increases, by a factor of 5 at 30 T and p:^20 at (the 
admittedly impractical) 100 T. 

3. G4beamline 

We next performed a similar study with the G4beamline simulation package. The "knob" 
comparable to absorber subdivision is the choice of the maxStep parameter, which is 100 mm 
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Figure 2. Root-mean-square value of x after 30 cm LH2 absorber vs number of slices, in 
Gaussian model, for longitudinal B fields indicated. 
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Figure 3. Distributions in x after 30 cm LH2 absorber with S = 30T, in G4beamline 
simulations, with (left) the default value maxStep= 100 mm and (right) maxStep= 1mm. 

by default. It represents the largest step along the particle trajectory that Geant 4 is 
permitted to take. Figure [3] compares the r.m.s. x distributions in a 30 T longitudinal field 
after traversal of a 30 cm LH2 absorber by 200MeV/c muons with maxStep = 100 or 1 mm. 
The default value of maxStep leads to an overestimate of the r.m.s. spread by about 50%. 

4. ICOOL 

Finally we studied the same situation using ICOOL. The relevant knob here is the choice of 
the ZSTEP parameter, which we took as 10~^ or 10~^m. In Fig. [i] we see that the spatial 
spread is independent of the value of ZSTEP. 

5. Discussion and Conclusions 

We have verified Lebrun's conclusion that high magnetic fields suppress some of the effects 
of multiple scattering and can thus produce smaller beam emittances in ionization cooling 
than one might naively expect. Because of this effect, G4beamline cooling results at high 
fields are not reliable unless the maxStep parameter is chosen carefully. We have discussed 
this with Yonehara, who used G4beamline to simulate helical cooling channels at fields up 
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Figure 4. Distributions in x after 30cm LH2 absorber with S = SOT, in ICOOL 
simulations, with (left) ZSTEP = 10"^ m and (right) ZSTEP = 10"^ m. 



to ^15T [7\. We found that he performed a maxStep sensitivity study and so was not 
misled by this effect. In contrast, ICOOL uses a variable-stepping algorithm by default, 
and so is largely insensitive to the user's choice of ZSTEP. Thus Palmer's simulations of 
30-40 T "final-cooling" channels [3J are also insensitive to the effect. Modification of the 
Geant 4 variable-stepping agorithm to take this physics into account is under study. Until 
it is successfully implemented, G4beamline users are warned to include a study of stability 
with respect to variation of maxStep if they are using fields of ^ 5 T and above. 

It is worth noting that the three simulations disagree at the ±10% level on the amount 
of scattering in 30cm of LH2, with (for example) Xr.m.s. — 0.38, 0.32, and 0.35mm (with 
neghgible statistical uncertainty) for S = 30T in the Gaussian, G4beamline, and ICOOL 
cases. This is not surprising, since the Bethe calculation [2J and the PDG formula (Eq. [T]) 
are based on approximate (Thomas-Fermi) atomic potentials, whose use is unnecessary 
for hydrogen [8j. G4beamline and ICOOL use more recent scattering models based on 
exact hydrogen wavefunctions [9j. The residual 10% disagreement between G4beamline 
and ICOOL was not investigated. 
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